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I  INTRODUCTION 


With  the  advent  of  high  speed  computers,  new  approaches  to 
the  solution  of  engineering  problems  have  arisen.  One  such  field 
is  Computational  Fluid  Dynamics,  in  which  numerical  techniques 
such  as  finite  differences  are  used  to  numerically  integrate 
partial  differential  equations  governing  the  physical  phenomena. 
One  problem  with  this  approach  is  that  on  the  boundary  of  the 
numerical  grid,  artificial  boundary  conditions  must  be  imposed, 
usually,  these  conditions  are  based  on  a  linearization  around  the 
flow  at  infinity,  and  are  thus  non-physical.  in  this  paper,  we 
perform  a  series  of  numerical  experiments,  in  which  we  evaluate 
various  boundary  conditions  of  this  type,  and  investigate  whether 
the  linear  model  is  an  accurate  representation  of  the  fully  non¬ 
linear  equations. 

The  flow  of  fluid  around  obstacles  in  two  dimensions  is 
described  by  the  compressible  Navier-stokes  equations 
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e  =  c^T  +  (u^  +  v2)/2 


Here,  the  four  variables  p,u,v,e  represent  the  physical  ^ 
quantities  of  density,  x-  and  y-  components  of  velocity,  and 
internal  energy. 

This  nonlinear  system  of  mixed  parabolic-hyperbolic  type  in 
two  space  dimensions  and  time,  with  four  independent  variables 
must  be  solved  in  an  exterior  region  in  .  The  geometry  will 
depend  on  the  particular  physical  situation  that  one  is 
attempting  to  model. 

the  situation  we  shall  be  interested  in  occurs  in  modelling 


Con  y 


flight  conditions,  in  which  the  conditions  at  infinity  are  pre¬ 
scribed  with  a  large  u-velocity  and  v-velocity  zero  and  pre- 
scibed  p  and  e  .  Fluid  flows  around  an  obstacle  in  x-y  space. 
Usually,  this  equation  is  solved  numerically  using  a  finite 
difference  scheme  of  the  Lax-Wendroff  type,  such  as  the 
MacCormack  ADE  method.  Since  these  calculations  can  only  be  made 
on  a  finite  grid  of  points  in  xy  space,  an  artificial  far-field 
boundary  is  created.  This  boundary  ought  be  sufficiently  far 
away  from  the  object  around  which  the  fluid  is  flowing  so  that 
local  phenomena  are  not  neglected  by  omitting  part  of  the  region 
of  fluid  flow,  on  the  other  hand,  the  farther  away  the  region 
is,  the  more  grid-points  need  to  be  included  a^fid  thus  the  more 
expensive  and  time-consuming  the  computations? become. 

One  then  has  the  problem  of  deciding  what  effect  this  new 
boundary  has  on  the  solution  of  the  problem.  Because  of  the 
viscosity  terms  in  (1)  and  the  additional  artificial  viscosi  ty 
introduced  by  the  finite  difference  schemes,  some  boundary  condi¬ 
tions  must  be  imposed. 

As  we  shall  show  in  this  report  on  numerical  experiments, 
considerable  care  must  be  exercised  in  the  choice  of  the  boundary 
conditions.  If  one  is  interested  in  steady  state  flow,  then  one 
starts  off  with  an  initial  approximation,  and  hopes  that  the 
errors  in  the  numerical  solution  propagate  out  of  the  region  as 
transitory  disturbances  in  the  physical  variables.  One  then 
expects  to  converge  to  the  steady  state  flow. 
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We  have  previously  shown  that  the  incorrect  choice  of 
boundary  conditions  can  give  rise  to  some  of  the  following 
phenomena:  i)  reflecting  boundary  conditions,  in  which  the 

disturbances  in  the  physical  variables  represented  by  the 
difference  between  the  steady  state  flow  and  the  initial 
conditions  are  not  allowed  to  exit  through  the  far-field  boundary 
but  instead  continue  echoing  within  the  grid,  and  giving  rise  to 
spurious  oscillatory  solutions;  ii)  under-specified  boundary 
conditions,  in  which  large  errors  are  introduced  before 

convergence  takes  place. 

# 

In  this  paper,  we  first  discuss  the  theory  for  simple  linear 
hyperbolic  systems  in  one  and  two  space  dimension.  We  then 
analyze  several  computational  experiments  in  the  light  of  this 
theory  in  two  space  dimensions. 


II  A  REVIEW  OF  THE  LINEAR  CASE 

The  method  used  in  the  calculations  which  are  the  subject  of 
this  paper  is  the  MacCormack  alternating  direction  explicit 
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scheme.  [6J  [15].  This  is  a  multistep  efficient  scheme  which 

reduces  in  the  linear  case  to  the  Lax-Wendroff  scheme.  For  a 
diagonal  N  x  N  matrix  A,  this  scheme  approximates  the  equation 
Ut  +  AUjj  =  0  by 


(2) 
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As  usual,  represents  the  space  step  and  n 

represents  the  time  step.  If  we  are  considering  the  equation 

f 

Ut  +  AUjj  =  0 

on  the  region  {(x,t),  0<x<l,  t>0}  then  the  'analytic  solution  is 
determined  by  the  initial  conditions  and  boundary  conditions  at 
X  =  0  and  x  =  1,  if  the  first  k  eigenvalues  are  positive  and  the 

remaining  N  -  k  are  negative  then  the  quantities  u^ must  be 

prescribed  at  x  =  0  and  - u^  must  be  prescribed  at  x  =  1, 

Thus  if  Wj  =  (Uj^,  \i2, - u,^,  0  0 - 0)  and  Wjj  =  (0,  0, - 0, 

- then  for  well-posedness,  the  boundary  conditions  must 

be 

Wj  =  f(t)  +  BqWjj  at  X  =  0 

‘  Wii  =  g(t)  +  BjWj  at  x  =  1.- 


This  gives  a  total  of  N  boundary  conditions.  If  either 
kx(n-k)  matrix  and  the  (n-k)xk  matrix  Bq  are  non-zero,  then 
the  boundaries  are  reflecting ,  that  is  a  wave  in  Wj  travelling 
left  to  right  will  be  reflected  as  a  wave  in  Wjj  running  right  to 
left. 

If  a  boundary  is  supposed  to  be  non-physical,  then  it  should 
not  be  reflecting,  since  the  reflections  would  depend  on  the 
location  of  the  artificial  or  numerical  boundary. 

NOW  let  us  consider  the  difference  scheme  of  (2).  If  the 
grid  points  are  given  by  {Xj}Q  with  x  =  0  and  Xj  =  1,  then  it  is 
clear  from  (2)  that  2N  boundary  conditions  are  required.  Thus  we 
must  prescribe  boundary  conditions  in  such  a  way  as  to  least 
affect  the  closeness  of  the  numerical  to  the  analytic  solution. 

This  situation  has  been  the  subject  of  several  papers.  In 

[3] ,  Gustaffson  and  Kreiss  point  out  the  danger  of  over- 

> 

specification.  By  this  is  meant  that  all  uj^  are  specified  at 
both  endpoints.  This  might  be  tempting  because  one  might  argue 

that  if  the  boundary  onditions  u^^ - Uj^  are  given  constants  at 

X  =  0,  then  eventually  these  same  values  will  be  assumed  by  these 
variables  at  x  =  1.  However,  in  [4]  it  is  pointed  out  that 
convergence  may  or  may  not  occur,  depending  on  whether  the  number 
of  grid  points  is  odd  or  even. 

A  method  which  works  well,  as  pointed  out  in  [3] ,  is  to 
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(0,t)  =  0 


k+l  <  i  <  n. 


or  in  the  numerical  scheme 


'^i » J  ^i  /  J— 1 


1  <  i  <  k. 


u. ,0  =  u .  ,  k+l  <  i  <  n. 

This  introduces  small  errors  at  the  outflow  but  these  errors 
do  not  propagate  upstream.  This  is  proved  analytically  by  Parter 

[10].  I 

In  [12] ,  many  different  numerical  boundary  conditions  are 
given.  The  conclusion  is  that  upwind  differencing  at  the  point 
of  outflow 


is  most  accurate,  although  it  converges  with  the  same  speed  as 
the  previously  discussed  u^^  =  0. 

The  most  serious  error  which  one  could  make  would  be  to 
prescribe  conditions  at  the  wrong  end.  In  other  words,  since  Uj^ 
is  right  running,  this  would  involve  prescribing  Uj^  at  x  =  1  and 
imposing  Uj^jj  =  1  at  x  =  0.  This  would  result  in  convergence  to  a 


steady  state  which  depends  on  the  initial  conditions 


III  THE  NAVIER-STOKES  EQUATIONS  AND  CHARACTERISTIC  VARIABLES: 

We  now  begin  our  discussion  of  the  equations  of  gas 
dynamics.  We  will  neglect  viscosity  for  the  purposes  of  this 
analysis.  We  will  assume  that  the  flow  is  one-dimensional  and 
subsonic  and  that  the  deviations  from  free-stream  solutions  are 
small.  This  will  allow  us  to  neglect  second  order  terms. 

There  are  many  forms  of  this  equation,  but  the  one  most 
suitable  for  the  present  discussion  is 


or  in  terms  of  physical  variable 


where 


A  =  am 

and 


Here  we  make  the  key  assumption  that  deviations  from  the 
free  stream  are  going  to  be  sufficiently  small  that  we  can  treat 
the  entries  in  the  matrix  A  as  being  approximately  constant  (at 

least  locally).  Denote  these  frozen  variables  by  0-subscript. 

0 

We  then  make  the  substitution 


(4) 


-  W  \ 
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and  when  this  is  substituted  into  (  )  we  obtain 


(5) 
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Notice  now  how  this  breaks  down  into  two  separate  cases.  On  the 
one  hand,  if  flow  is  supersonic  then  all  wave  motion  is  in  the 
left  to  right  direction.  In  this  case  all  analytic  boundary 
conditions  should  be  prescribed  at  the  left  hand  side  and  only 
numerical  boundary  conditions  prescribed  at  the  right  hand  side. 
Since  the  substitution  (4)  is  equivalent  to 

p  =  +  ( Pq/2Cq) (K2+K3) 

u  =  1/2(K2-K3) 

it  follows  thfet  prescribing  all  physical  variables  at  the  inflow 
and  prescribing  3p/3x  =  3u/3x  =  3p/3x  =  0  at  the  outflow  is 
legitimate  in  terms  of  analytical  and  numerical  requirements  in 
the  supersonic  case. 

However,  we  must  now  consider  the  case  of  subsonic  flow.  in 
this  case  the  situation  is  completely  different.  Here,  two  of 
the  variables  and  W2  go  left  to  right  with  velocities  Ug  and 
Ug  +  Cg  respectively,  whereas  one  of  the  variables  runs  right  to 
left  with  velocity  Cg-Ug .  The  variables  W2  and  W3  have  no  clear 
physical  significance,  yet  it  is  only  in  considering  these  variables 
that  the  full  wave  structure  of  the  equations  (5)  or  (6)  can  be 
understood.  Thus,  one  would  be  led  to  predict,  for  small 

deviations  from  free  stream  conditions,  that  the  best  boundary 

/ 

conditions  would  be,  for  an  interval  (0,L) 


(0,t)  = 


{L,t)  =  0 


W2  (0,t)  =  K2 


(L,t)  =  0 


(0,t)  =  0 


(L,t)  =  K3 


Note  the  curious  aspect  of  these  boundary  conditions.  In 
order  to  prescribe  the  numerical  values  and  K2»  we  need  to 
know  accurately  all  three  physical  variables  at  some  distance  to 
the  left.  However,  only  the  two  combinations  and  K2  are  pre¬ 
scribed.  This  can  be  summarized  by  saying  that  while  we  have 
used  all  three  pieces  of  information  upstream,  we  have  done  so  in 
such  a  way  that  one  degree  of  freedom  remains^  thus  allowing  the 
waves  in  W3  to  exit  without  reflections. 

On  the  basis  of  the  linearized  model,  various  other  combi¬ 
nations  would  be  well-posed.  For  example,  \t  is  possible  to 
prescribe  K3  in  terms  of  either  or  K2  at  the  outflow  x  =  L. 
Thus  at  the  outflow  one  may  prescribe 


W2  (L,t)  =  F3  (t)  +  c^W^  (L,t)  +  C2W2  (L,t) 


For  example  if  =  0,  C2  =  1»  then  this  amounts  to  putting 


u  (L,t)  =  1/2  F: 


(i.e.',  we  prescribe  velocity  at  the  outflow) 


Alternatively,  we  might  take  Cj  =  0,  C2  =  -1  and  we  would  get 

(8)  p  (L,t)  =  ((PoCo)/2)F3 

(i.e.  we  prescribe  pressure  at  the  outflow). 

Many  other  combinations  are  possible,  but  as  remarked  in  section 
II,  all  these  will  cause  errors  in  the  initial  data  to  be 
reflected  back  into  the  medium  as  waves  running  from  right  to 
left.  For  example,  we  would  predict  that  an  error  in  W3  would  be 
reflected  back  as  an  error  in  if  we  use  boundary  condition 

(8) .  As  we  shall  see,  this  is  exactly  what  happens. 

At  the  inflow  end,  we  may  prescribe  and  W2  in  terms  of 

W3.  Thus  the  following  boundary  conditions  are  well  posed; 

# 

r 

(9) (a)  Wj  (0,t)  =  F^  +  CJW3  (0,t) 

(b)  W2  (0,t)  =  F2  +  C2W3  (0,t) 

/ 

For  example,  choosing  C2  =  +1  in  (9b)  corresponds  to 

u(0,t)  =  (1/2)F2 

(i.e.  prescribing  u  at  the  inflow)  and  C2  =  -1  corresponds  to 

p(0,t)  =  (PoCo/2)F2 

(i.e.  prescribing  p) .  One  can  prescribe  the  combination  (u,p)  by 
first  choosing  C2  =  1  (thereby  prescribing  u)  and  then  choosing 
Cl  *  Pft/Cftf  thereby  prescribing  p  in  terms  of  a  given  F,  and  a 


o 


VI  TWO  DIMENSIONAL  RESULTS 

Having  now  understood  the  phenomena  which  can  occur  when 
calculations  are  made  with  one  space  dimension,  we  now  consider 
the  more  complicated  situation  of  two  space  dimensions. 

In  this  situation,  we  shall  solve  the  Navier  Stokes  equation 
(1),  again  by  the  standard  MacCorroack  A.D.E.  scheme  on  a  20  x  20 
grid,  which  is  a  very  simplified  model  of  a  wind  tunnel.  We 
shall  continue  to  consider  flow  close  to  the  free  stream  flow 
(M^  =  0.5)  previously  studied  in  the  one  dimensional  case.  The 
physical  values  are  given  earlier.  Figure  8  shows  the  geometry 
of  the  situation.  The  fluid  is  flowing  in  at  the  top  right 
corner  of  the  grid  and  flowing  out  at  the  bottom  left. 

We  have  three  distinct  types  of  boundaries  to  consider.  We 

# 

have  the  inflow  and  outflow  boundaries  (as  bfefore)  and  in 
addition,  two  sidewall  boundaries,  where  the  fluid  is  flowing 
parallel  to  the  boundary.  ^ 

We  shall  continue  to  impose  one  dimensional  boundary 
conditions  of  the  type  given  in  the  first  section  on  the  inflow 
and  outflow,  along  with  the  additional  condition  v  =  0.  This 
says  the  fluid  flow  is  one  dimensional  at  the  inflow  and  outflow, 
and  could  be  physically  reasonable. 

When  we  come  to  the  sidewall  conditions  we  must  undertake 
another  one-dimensional  analysis.  Thus,  we  assume  all  variables 
are  constant  in  the  x-direction  and  variation  only  *akes  place  in 

the  y-direction. 
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This  leads  to  the  set  of  equations 

(6)  Vt  +  AVy  =  0 

where 
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Freezing  the  coefficients  of  A  in  (6),  we  substitute 


(7) 


T,  =  p  +  ^ 


=  V  +  - ^ 

2  PnC 
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3  PnC 
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=  u 


Equation  (6)  is  then  transformed  to 


Thus,  if  we  only  consider  deviations  transverse  to  the  free 
stream  flow,  we  have  on  the  basis  of  the  linearized  model,  four 
non-physical  variable  ,  T2 »  T3,  T4 ,  two  of  which  are  u  and 
entropy,  Tj  move  with  zero  velocity  in  the  y  direction,  one  of 
which,  T2,  moves  with  speed  Cq  in  the  positive  x  direction,  and 

one  of  which  moves  with  speed  Cq  in  the  negative  x  direction. 

# 

Thus  one  dimensional  theory . predicts  that  at , the  sidewall,  we 

# 

should  impose 


y  =  0 

Tiy  =  const 

T3y  0 
T4y  =  0 


This,  together  with  the  inflow  and  outflow  conditions  gives  the 
following  set  of  boundary  conditions 


Inflow 


Vj  =  V, 


Pi  ”  4  'P2  *  <'0=0*'‘3'“2” 


“1  =  I  '^3  -  7^  *  “2> 


Pi  =  4  ^  ^  '7^  ^  S  -  ''2> 


where  the  free  stream  values  are  specified  in  the  characteristic 
variables  kj  =  e,  -  P./Cq^  kj  =  U,  +  P./PoC^  and  where  the  zero 
subscript  refers  to  frozen  variables 


Outflow 


''n  =  '^n-l 


“n  ■  5  ‘“n-1  *  Pn-l/'Po=o’  '  "4 


Po‘=0 
Pn  ”  — 2“ 


"‘4  *  "n-1  Pn-l/'Po"o” 


n  =  O 


^  [k..  +  u.  ,  + 


Top  wall 


»0=0 


£i-i 


Pj  "  7^  ■  "s' 


Vj  =  l/Zkj  +  l/2(v^_j  +  Pj.j/Po^O^ 


P. 


P-5  = 


(V, 


+  liii- 


j  2Cq  j-1  PqCq 


•^3^  Pj-1  "  Pj- 


Bottom  wall 


ui  =  li«o 


Vj  =  1/2[v2  +  kj  -  P2/^P0^0^^ 


Pi  “  — 2“  ^“'^2  '''  ^2  ^ 


“l  *  ‘’2  ‘V^^o’  '■''2  ''2  ■  V‘'o‘'o' 


where  the  free  stream  values  are  used  to  specify  the 


characteristic  combinations 


The  purpose  of  the  present  grant  was  two-fold:  first  to 
expand  the  previous  work  to  situations  where  the  flow  was  at  an 
angle  a  to  the  artificial  boundary,  but  is  primarily  steady- 
state,  and  second,  to  see  how  the  artificial  boundary  conditions 
apply  in  period  situations  where  convergence  to  free  stream  is 
not  expected. 

The  first  objective  is  clearly  desirable,  simply  to  save  in 
additional  computing  time.  The  requirement  that  the  grid  not  be 
rectangular  with  boundaries  either  parallel  to  or  perpendicular 
to  the  free  stream  at  infinity  clearly  can  be  used  to  remove  a 
large  number  of  points  from  the  grid.  The  entire  boundary  can  be 
chosen  approximately  equidistant  from  the  object  whose  flow  char¬ 
acteristics  are  being  studied. 

The  second  objective  was  chosen  because  of  its  importance  in 
modelling  flow  in  the  compressors  of  turbines.  If  one  is 
attempting  to  design  turbine  compressors,  one  must  have  knowledge 
of  how  the  flow  in  the  area  between  the  stators  and  the  rotors  is 
behaving.  This  is  almost  impossible  to  measure. 

On  the  other  hand,  to  accurately  model  the  entire  inside  of 
the  compressor  requires  a  computational  capability  far  exceeding 
the  largest  of  today's  computers.  Thus  one  models  only  the  flow 
past  two  blades  of  the  rotor,  and  uses  periodicity  to  extend  the 
calculations  to  the  entire  rotor.  The  flow  coming  in  past  the 


o 


stator  is  regarded  as  a  periodic  inflow  term. 

The  following  pages  summarize  the  results  of  the  research  on 
these  two  projects. 

Section  I  Flow  at  an  angle  of  incidence  a. 

The  previously  outlined  boundary  conditions  were  used. 
Computations  were  done  on  the  basis  of  the  inflow  or  outflow 
depending  on  the  flow  perpendicular  to  the  surface.  For  computa¬ 
tional  simplicity,  we  chose  a  rectangular  grid,  but  choose 
u  a  V  *  300  This  results  in  a  45®  angle  of  incidence  or 
exit.  We  impose  an  initial  10%  disturbance  in  W2,  the  right 
running  characteristic  variable  and  observe  what  happens  as  it 
exits.  Recall  that  this  disturbance  runs  at  an  angle  of  45®  to 
the  flow.  The  following  plots  for  W2  at  various  iteration 
numbers  show  it  exiting  without  any  surprising  phenomena,  similar 
to  the  one-dimensional  case.  By  the  sixty-fifth  iteration,  the 
flow  has  almost  converged  to  steady  state.  There  appear  to  be 
some  small  reflections  (as  shown  in  the  last  plot)  in  W3,  but 
nothing  to  worry  about. 

Thus  we  can  conclude  that  the  inflow-outflow  boundary 
conditions  work  reasonably  well.  Subject  to  a  rotation  of 
coordinates,  they  can  work  on  any  geometry  -  the  only  decision  is 
whether  we  are  at  an  inflow  or  outflow  point. 


II.  The  case  of  periodic  BCs 

We  now  consider  some  experiments  on  the  case  of  periodic 
boundary  conditions.  We  study  flow  parallel  to  the  rectangular 
grid,  perpendicular  to  the  inflow  and  parallel  to  the  side-walls. 
The  most  obvious  boundary  condition  is  pure  periodic. 

=  S2  foe  all  physical  variables  S.  This,  it  was 
feared  might  lead  to  over  specification,  but  these  fear  appeared 
groundless.  A  sight-traveling  wave  would  exit  to  the  right,  and 
reappear  at  the  left  in  a  predictable  and  correct  manner. 

The  next  objective  was  to  see  if  this  could  be  done  by  only 
using  periodicity  in  the  two  sideways  traveling  periodic  waves. 
The  reason  for  this  was  to  minimize  storage  problems  for  the  more 
complicated  turbine  problem  already  discussed. 

Curiously,  this  method  failed.  Furthermore,  the  failure 
seems  intrinsic  to  the  whole  idea,  as  we  shall  demonstrate. 

First  the  W2  wave  appears  in  iteration  5  and  proceeds  to 
exit.  Periodicity  requires  that  it  reappear,  and  this  can  be 
seen  in  iteration  20.  However,  already,  a  sharp  spike 
(reminiscent  of  over  specification  is  appearing  in  W3  by 
iteration  20.  As  far  as  iteration  40  W2  continues  to  progress 
from  right  to  left,  but  at  iteration  55,  sharp  spikes  begin  to 
develop  in  W2  and  W3 .  This  seems  to  suggest  overspecification  at 
the  corners,  but  as  far  as  we  know,  this  is  not  the  problem. 
Shortly  afterwards  the  program  crashes.  The  relevant  portion  of 
the  program  for  the  boundary  condition  implementation  is 


included. 


ITERATm  NUMBER 


ITERATION 


ITERATION  NUMBER 


ITERATION  NUMBER 


ITERATI 


o 

II 


pc'. 

W 

pq 


o 


ITERATION  NUMBER 


ITERATION  'NUMBER 


«« 

- 

M 

c 

44 

N 

^rtj 

4  -J 

n 

«• 

lA. 

V  W 

4 

w 

w 

•4M 

Z 

ft  W 

•7 

Al 

c 

« 

44^ 

f  o 

4 

M 

3 

V4 

4 

^  o 

M 

1 

« 

•  « 

«  ^ 

4^  X 

W 

4  4 

«« 

♦  fti-> 

-o 

4  ft 

44 

V  « 

« 

U. 

ruz 

1  N 

o 

w  «4 

w 

z 

O.— 

w  4«4 

4  ^ 

3 

• 

O  4 

w 

44 

b.  Ow 

N  O 

C  44 

O 

4 

Z  c 

a. 

V 

z  zo. 

lA.  m 

•J 

Z 

w 

c 

ft  ^ 

:> 

m 

•-*  <k5 

zo 

V  1 

ft 

# 

W 

\  “7 

o 

u 

O  NO 

«Z 

4^  4 

3 

4 

-^x 

z 

m  ^x 

044 

o  c 

♦ 

O 

4<  V 

mac 

M 

U  ^OE 

XO 

o  ^ 

X 

^  V 

4  «4 

m  w 

o 

Z  ^  %W 

ftZ 

a; 

ft 

c  o 

X 

•-•  #  1  ♦ 

'-ft 

ft  w 

«« 

'O 

•n  X 

1^  <• 

tt 

4  w 

N  O 

«» 

4 

niz 

1  ft 

z  ai 

w 

Z  OIOI 

N 

3 

'ft 

4  N 

« 

N 

oc  «  «•  * 

cu« 

«4  O 

44 

cu 

V  s 

V  — 

OM 

a  «* 

4  Z 

4 

M 

•4  44 

-» 

N  3-^ 

M  4 

1  ft 

1 

M 

04« 

3  * 

N  *4 

1a.C 

4 

44 

X  • 

O  C 

^  « 

01 

-1  X  -* 

Z.J 

C  4- 

3 

44 

KC 

X  -J 

*  O'-^ 

4-1^ 

-J  ▼ 

V 

• 

n3 

ft  1 

Al  « 

w» 

a 

fU  -»  % 

o  w 

V  V 

w 

-7 

4  4 

*W 

X 

O. 

^  1 

4^0 

w  O 

o 

44  4 

.*•  V 

W*- 

o 

w 

&  rj  w 

N  O 

a  w 

3 

V 

•V 

V  4- 

w 

O.  VO 

— *  w 

3  m 

O 

M 

M 

niw 

3 

a  3 

«« 

w  03 

4 

Q  u. 

X 

«* 

M 

3 

a  o 

»4 

1  — o 

Z  Z 

ft 

«* 

M 

o 

VO 

n  I 

1  X 

o 

n  MX 

C4-» 

ft  •*• 

w 

M 

X 

'■'X 

3  ft 

oc 

z 

«4 

V  OK 

M 

-J  4 

w  O 

M 

M 

ft 

3a 

M 

1  4 

o 

z 

«  Z  w 

M 

V7 

•*  •» 

4 

«• 

M 

O  4 

IP 

M 

w  ^ 

•  w 

o 

44 

4-  44  W 

M 

4-  4 

-4  O 

U 

«• 

M 

I** 

Z 

M 

M  W 

^  ♦ 

M  0  4 

M 

OC 

-4  Z 

ft 

M 

M 

ftO 

o 

M 

>*-  o 

> 

«• 

■••O 

M 

o-j-^ 

4^414 

N 

MOW 

1 

I  z 

w 

44  44 

a  or 

or 

•  ^ 

--^  Oft 

«• 

f  V-H'7-'0 

MUiM 

^4 

1- 

M 

a  N 

« 

OJ  — 

Z  N 

•» 

^  W  4 

3 

«t  UJ 

M 

44 

wo 

44 

^C'^3 

o 

«« 

w  • 

*•  4-*  44  <— 

M 

-40*7 -JOft 

O 

M  3  « 

4 

ZM 

o 

44.10'- 

4  3 

z 

« 

N  *7 

4^03 

«• 

•  O  4 

V  Z  4 

••  -4 

W  tf>  w 

"7 

l-tli. 

z 

•1Z  M 

■^O 

3 

M 

h.  • 

•XO 

-7XC 

w  44  A 

M 

M 

or 

o 

O 

C  'MO 

O 

zru 

4  4“»or  m 

•ft  ^o  o-> 

■7  • 

W  ft  «• 

V 

o 

-i 

.i  voz 

A 

•• 

M  M 

~im 

c 

C  N  V  O  ••  4 

«<U 

MOM 

o  o 

in 

1  — M^4 

Al««  • 

«• 

oo. 

•wf\|«-^ 

ft 

^  W 

XO  c 

•J  •* 

M  H 

aj 

zz 

> 

44 

-ooo 

C 

«« 

X  o 

an 

V  *40arz.j 

V  ** 

MOW 

»-*ft 

ft 

o 

VXZM 

a 

•• 

(k  X 

’-'OO  *-'49 

ft 

w  *3  «  •-•  V 

W  —4 

MIaJ 

I  ‘"‘ 

o  — 

ft 

wftP-fO 

w  —  ^ 

lii 

-'0:0  X  3  ni  • 

K 

O-iOlPO  w 

O  4 

M  ft  M 

ft  •"* 

X  N 

o 

V 

3  WOZ 

0.0.  — 

or 

c 

♦  \ 

ofto  ^aoi 

ap 

O  4X 

*ZO 

0*7 

M  ft  M 

ft  •*“ 

z 

ft 

oinx»4 

oo  • 

»- 

X 

w4 

Z  N  XO  0.  «• 

z 

zcaovo 

4 

MO 

M  4« 

c 

3 

o 

X  -ftO 

X  -  — 

Ml 

ft  ft  o  - 

3 

q:«j  • 

•  W  w 

•  -J 

MOM 

^  4 

N  -J 

o 

ftO'-X 

q:  •  « 

z 

o 

X  • 

O 

•  V 

'^  \  m 

V 

M 

M 

3*7 

#«*  -7 

M 

o 

•  •  N  ft 

•  -*  -/ 

3 

•  *» 

o 

^  ^  9H 

---'aA 

«• 

M 

^-JO  •> 

.J<\l  • 

o 

—  w 

O 

»/t 

-J*“  ^  •-•  “> UJ 

-40  4 

4*7  • 

•a 

M 

M 

*7  V  N  V 

V  'W 

an 

H- 

V*4  44CIA 

-.  .a 

£• 

O 

n-^fu 

•3*71  40 

10 

M 

M 

4  •O'* 

lA 

4  4  4^  • 

•  — 

3 

•  *  » 

H 

lO  4 

4C  ■ 

♦x 

o 

M 

M 

44  44  ft  O 

ft 

44^  J1» 

-«  «T  UJ 

4  ^  4  4  «  ^ 

m 

-^oan 

w  •  •  w  T  an 

•  o  •  an 

UJ 

w 

•  11  •  • 

a  3 

4* 

•  *3 

-.40. 

«f 

-JOT  W  V7-* 

V  -  3 

44 

31V  •  ft  3 

V  •  ~  3 

M 

V  •  •W'^ 

r 

">  —  ‘3 

.  3-->  -o 

M 

V  4 

w  .  %<* 

4^  N  z 

r  -“wr 

''  O’  Z 

aP 

W 

Ww-J 

-  a 

f* 

-  x  a  oo  -a  X 

.400  O  3-’ 

O  *4  44 

o 

(-•\au>i9*^ 

afiV*-  ••« 

r 

M 

44  S4  W01 

O  Mi  ' 

bf 

♦no  Or3ZX  '-«MiO. 

«» 

O  V3 

30  V44 

W-4>- 

»- 

tM  —  .Jl- 

ft 

M 

4*  33  4  4 

C  C-  '  T 

O 

«» 

o  - 

oo  o - o  - 

*• 

0X0 

OOlZ 

Z  V(UZ 

—  iZ 

ft 

OO  *»  W 

X  X  MC* 

ox  1 

X  -  XCl  X 

«• 

X  -^1140*7100 

o 

OOO w . o 

orxuo 

•- 

OIX-4W 

O  O. 

#/ 

<* 

CiO 

tf  ft  o  o 

ft  ft  or  o 

ft  o 

o 

oooo  o 

O  3  30 

Oft  ft  'ft 

•• 

o 

m* 

^4  «4 

«* 

—4 

—4  —4 

— 

44 

44 

o  ao 

aa 

rj 

rt 

v4  44 

<u 

V-  Ow  oooooo  o  ooooo 


^  —  r:'*  ^  j  fc'f  !  c  c  '••Aj^  '*4^TU’''rf-oo»^— cdo)9  ^ni 

uuu;;t/u«.b'*t«/  i‘L  t0  i  ‘C  -4.  i-t  r•f•f•r•r•r*f■r'rwoL.O)0}<3B^7(9(r«*voo•«f*•k^9>(^c^c^c^r^c^o90A«»SA«^o<^«-•-••« 


!  Pi >  .11’  •  ;nf  ' ) 

1,1-1  4:  CU  -  .  •  .<1.1  .  .PM.iii »  .Jiff,  t  JL  >-Pl*r,  JL«  )  >-'(CINF*»2  > 

'••••  -H  E  ■  .  I  O' k  ..'I '«C^' IBHOUIK.  JL.J  )t»2)*tRH0U(K,JL.  n* 


3 

4 


<U 


3 
o 
X 
^  a 


o 

r 


a; 

w 


PH  ^  o- 


O  HP  , 


fU 

w 


(/I 

z 

o 


•-K 

« 

M 

r» 

E 

• 

V 

PHfU 

1*4 

▼ 

pJ 

^  % 

« 

X 

Z 

w 

W 

* 

4 

•■H  PH 

« 

w 

“7 

fua 

0^ 

1 

0 

-P  “7 

M 

pH 

X 

pj 

^  « 

u. 

«  « 

•  Ui 

•<4 

Oe 

V 

h>*4 

Z 

<u-> 

niz 

• 

w 

am 

A 

-7 

3 

•-7  * 

0 

u. 

0*4 

N  0 

« 

0 

M 

z 

Xw 

idk  m 

C 

» 

X 

p^-jcr 

id. 

4X3 

XU 

.J 

•7 

CX 

*  *> 

z. 

0 

N  0 

>^z 

W 

% 

w 

9>  « 

«• 

^  X 

o*^ 

c 

4 

ruou 

0 

u. 

^OL 

0 

■ij 

PH 

H  *Z 

X 

z 

fk  X 

X 

ru 

ax»-* 

a 

•-» 

“7  ♦ 

K 

M 

ruoQ. 

0 

*  w 

S. 

3 

•• 

^  *  H 

s 

X 

oin) 

S 

0 

-• 

a  K 

(X 

w  M 

ru  '•H 

•>• 

X 

UtOCL 

-> 

3M 

*»  -7 

K 

• 

*  »  % 

N 

0^ 

M  » 

-> 

w 

-> 

<t  3  > 

ru 

X 

Ul  Z 

4 

* 

C»~  X 

*-•  V 
Ohp 
wa 
S  w 


C 


o 


a  u.  c 

0£  V) 


3 

•  W  3 

z 

-  -  «3 

‘^4  Oft 

X 

0 

ru  H-o 

0 

^1  wO 

x  z 

ft 

X 

*  IL  Z 

0 

u.  0  «  I 

»«<«  ft  M 

w 

ft 

Z  V  ftft 

c  •  W  U 

> 

•-•  ft  Zw 

•• 

•J*7  •• 

4 

«• 

ft 

0  HP 

w 

w  •  ft 

0 

«• 

m 

« 

••  -H  ^  ^ 

w 

wC  z 

ft 

0 

•  •*  U 

M 

aj.-  %— 

X 

« 

z 

« 

ru--*  •  —  u.  ft 

W 

♦  ^  •-••i  —  0 

— • 

<-r  ft  S 

3 

«• 

HP  -z,  N 

M 

—  W  *  HU.  X 

3 

0 

«« 

S  ^ 

«« 

-•0-7-JZ  ft 

0  “* 

— 

1*.  -  ♦  3 

m 

«• 

U*7'^-»^03 

•» 

•  X  hW  ^  4 

«*  <9 

<UWU.Z'- 

M 

Z  •  hXO 

M 

->tt  t  H-O- 

PH  • 

•  —  z 

c 

«» 

— •'V  •rv“>«  •• 

M 

•  N  -lO  ••  “1 

-7fU 

ru 

_J 

W'T/h- 

ft 

«• 

0-*p-7w  «<•>  PH 

«* 

E  ^  X  ft  • 

%«• 

't 

UI 

«* 

10  •OfU'-'-l  — 

.iJ^^wftZ  c 

P.J 

• 

3 

ft 

ft  I  -•  I  —  ^ 

V  •  3  •  ►-•3 

V 

p« 

3  «•  u 

#- 

.  wft  wft:>  . 

z 

w  H^OI/IOV 

H  H 

» 

z 

'jM.X 

VI 

c 

♦  NON  oru '-' ru 

ft 

0  *1  •  X  H- 

ft  '7 

pJ 

X 

0  z  oa  H 

z 

ft 

UJ 

XZftOftft 

w  • 

f- 

to 

I  *-,x  -  -x 

3 

UI 

V '^fk  — ftft  w- 

ft 

ft  ^  •  •  W  4P, 

•  •-> 

»- 

4.1  ■•  T 

a  Oft  «  H 

0 

ft 

—tft  *• 

«  ^  P«  .-4  S  M 

K  hX 

a.  u.  X 
-ZH 


ObL  H 


3  O*^ 
OZ  • 
•ftC 
«  •^i 
t-  M.  ^ 


. 

, 

.p 

l/<ft 

-  --  .w 

V* 

z 

P 

« 

p 

—  -o 

1^ 

a 

— •  p*  rtj  « lU  •  •  •  H 

3 

z 

■ 

- 

.-.  I  UJ 

3 

•  OH 

0 

% 

• 

^■1 

-•  pft  3 

r 

•7*7  3  -3  H  • 

• 

-  , 

* 

-  ru 

X 

•  3  -^0  -^0 

« 

p 

ul 

•• 

*3poThI4.->w1 

“1 

” 

p 

.  ud  ■  h- 

u* 

“•  H  ft  -J  ft  .•»  H|*>  ft 

«» 

>P/ 

.P 

ft 

^  pr 

o«:  -  r 

»- 

• 

ft 

\ 

p 

<1^ 

2 

1 

• 

2  2^0 

•p 

OIlX4X-Ix 

ft 

u 

3 

ft  u:  4  *. 

a 

•» 

oa  ft  ft  ft  a 

ft 

-• 

— 

•- 

4 

— •  p-  •- 

ft 

•<3  • 
•o-#“ 
-»X  * 
-a  p 
..i  ♦  V  ^ 
V  T  - 


.  U>  ••  -p 


•o 

“I  I 
'  *0:0) 
'Jh  d 
'Vi  '  z 


iMi-^  r 
•33^ 

r<x«  • 

KUiUI  W 

»-  oa  * 

Ui  V 

z  u  • 

■-»  V  v»«U 


ozo 

X  w  I 

a.  ix 


x  ♦  - 
a.  a 


^  w  *  ^- 

'  I  *«  O 
tk.  v> 


ZiZ 

oor  • 

a  r  4  -- 
4-  X  »- 

u<  .r  j  •  c*  V 
o  uis/'w' 


a 

3 


2  % 


^•O 


H  -  p- •*!»•»  J  Jf  •  I*  «  J  |r«  ^  Ij.  ^4.f  .Mt»»  ru*"  ▼  4/<  sp  r*  fVC^  •!»  —  <U 

•••r-T**  1  t  I  j.«  u  ..  M-u  u  u  li  ii 

c  f  !  ;  .  .  I  I  f  .  i-f  f  f  I  f  rrr  1  fr-f  i-f.fHf-r-r- 

-  2.  Zf\  V-  \  i  ?,  -SSA  S  &  \  ‘34»<»‘%*'‘va» 

VVA/«va 


ITERATION  NUMBER 


ITERATION  NUMBER 


ITERATION  NUMBER 


ITERATION  NUMBER  =  20 


L 

m 

isSfiWf 

IM' 


•»  •  o 
^  M. 


ITERATION  NUMBER  =  25 


ITERATION  NUMB 


ITERATION  NUMBER  =  35 


z  '  o 


ITERATION  NUMBER 


ITERATION  NUMBER  =  55 


TTERATI 


Finally,  we  attempted  to  run  periodic  characteristic  B.C.s 
on  a  one-dimensional  basis.  We  felt  that  since  the  problem 
manifested  itself  at  the  corners,  if  we  removed  the  corners  the 
problem  should  go  away.  We  next  show  what  happened  in  this 
situation. 

First  we  see  the  one-dimensional  wave  in  W2  moving  to  the 
right,  exiting  and  reappearing  by  the  iteration  20.  However, 
sharp  spikes  appear  in  W^  by  this  time.  W2  appears  to  behave 
well  until  iteration  85,  but  either  reflection  or  some  non-linear 
instability  lies  caused  W3  to  grow  catastrophically.  Clearly 
this  represents  failure  of  the  boundary  conditions,  even  in  the 
one-dimensional  case. 

Various  minor  variations  were  tried,  but  all  produced  the 


same  basic  phenomena. 


Conclusion 

Non-reflecting  characteristic  boundary  conditions  work  well 
so  long  as  the  flow  is  basically  returning  to  steady  state.  This 
was  seen  in  the  first  series  of  numerical  experiments. 

However,  at  this  stage  of  development  they  cause  strange 
instabilities  which  prevent  efficient  modelling  of  periodic 
phenomena.  The  reasons  for  this  are  not  understood  at  the 


present , 


BIBLIOGRAPHY 


[1]  Engquist,  B.»  and  Majda,  A.,  "Absorbing  Boundary 

Conditions  for  the  Numerical  Simulation  of  Waves",  Math 
Comp.  vol.  31,  pp  629-651,  (1977). 

(21  Gottlieb,  D.,  and  Turkel,  E.,  "Boundary  Conditions  for 
Multi-Step  Finite  Difference  Methods  for  Time  Dependent 
Equations",  J.  Comp.  Phys.  Vol.  26,  pp  181-196  (1978). 

(31  Gustaffson,  B.,  and  Kriess,  H.O.,  "Boundary  Conditions  for 
Time  Dependent  problems  with  an  Artifical  Boundary", 

Jour.  Comp.  Phys.  Vol.  30,  pp  333-351  (1979). 

(4]  Gustaffson,  B.,  Kreiss,  H.O.,  and  Sundstrom,  A., 

"Stability  Theory  of  Difference  Approximations  for  Mixed 
Initial  Boundary  Value  problems  II",  Math.  Comp. ,  vol.  26, 
pp  649-686  ( 1972) . 

(51  Kreiss,  H.O.,  "Stability  Theory  for  Difference 

Approximations  of  Mixed  Initial  Boundary  Value  problems 
L"»  Math.  Comp.,  Vol.  22,  pp  703-714  (1968). 

(61  MacCormack,  R.W.,  "Numerical  Solutions  of  the  Interactions 
of  a  Shock  Wave  with  Laminar  Boundary  Layer",  Lecture 
Notes  in  physics,  Vol.  19,  Springer  Verlag  (1976). 

* 

t 

(71  McKenna  P.J.,  Graham,  J.E.,  and  Hankey,  W.L.,  The  role  of 
far-field  boundary  conditions  in  numerical  solutions  of 
the  Navier-Stokes  equations,  AFWAL  Tech.  Report  TR.  81. 

(81  Moretti,  G.,  Comment  on  "Stability  ^^spects  of  Divergent 

Subsonic  Flow",  AIAA  jour.,  Vol.  19,\No.  5.,  p  669  (1981). 

(91  Oliger,  j.,  and  Sundstrom,  A.,  Theoretical  and  practical 
aspects  of  some  initial  boundary  value  problems  in  fluid 
dynamics,  SIAM  J.  of  Appld.  Math.,  35  (1978)  pp.  419-446. 

(101  parter,  S.V.,  "Stability,  Convergence  and  Pseudo-Stability 
of  Finite  Difference  Equations  for  an  Over-Determined 
problems",  Numer.  Math,  vol.  4,  pp  277-292  (1962). 

(Ill  Rudy,  D.H.,  and  Strikwerda,  J.S.,  "Boundary  Conditions  for 
Subsonic  Compressible  Navier-Stokes  Equations",  Computers 
and  Fluids,  Vol.  9,  pp  327-338  (1981). 

I12j  Rudy,  D.H.,  and  Strikwerda,  J.C.,  "A  Non-reflecting 

Outflow  Boundary  Condition  for  Subsonic  Navier-Stokes 
Calculations",  jour.  Compu.  Phys.  vol.  36,  pp  55-70 
^  (1980). 


113)  Schlichting,  H./  "Boundary  Layer  Theory".  Seventh  Edition, 
McGraw  Hill  N.Y.  (1979). 

[14]  Shang,  J.S.,  Oscillatory  compressible  flow  round  a 
cylinder,  proc.  AIAA  20th  Aerospace  Meeting,  jan  11-14 
(1982).  Orlando,  Florida. 

[15]  Shang,  J.S.,  "Numerical  Simulation  of  Wing  Fuselage 
Interference",  AIAA  19th  Aerospace  Sciences  Meeting, 
AIAA-81-0048,  Jan  12-15,  1981. 

[16]  Steger,  J.L.,  Pulliam.  T.H.,  and  Chima,  R.V.,  "An 
Implicity  Finite  Difference  Code  for  Inviscid  and  viscous 
Cascade  Flow",  AIAA  Preprint-80-1427,  AIAAA  13th  Fluid  and 
Plama  Dynamics  Conference,  July  14-16,  1980. 

[17]  Yee,  H.C.,  "Numerical  Approximation  of  Boundary  Conditions 
with  Application  to  inviscid  Equation  of  Gas  Dynamics", 
NASA  Technical  Memorandum  81265,  March  1981. 

[18]  Yee,  H.C.,  R.M.,  Beam,  and  Warming,  R.F.,  "Stable  Boundary 
Approximations  for  a  Class  of  Implicit  Schemes  for  the 
One-Dimensional  Inviscid  Equation  of  Gas  Dynamics",  AIAA 
paper  81-1009,  presented  at  Palo  Alto,  CA,  June  1981. 


